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Spin relaxation in-plane anisotropy is predicted for het- 
erostructures based on zinc-blende semiconductors. It is 
shown that it manifests itself especially brightly if the two 
spin relaxation mechanisms (D'yakonov-PereP and Rashba) 
are comparable in efficiency. It is demonstrated that for the 
quantum well grown along the [001] direction, the main axes 
of spin relaxation rate tensor are [110] and [lIO]. 



I. INTRODUCTION 

Spin relaxation processes have significant effect in op- 
tical and kinetic properties of semiconductors. They play 
important role in optical orientation of electrons and nu- 
clei [yj and in anomalous magnetoresistance caused by 
weak localisation J2|. Both theoretical calculations and 
experimental data analysis have been carried out assum- 
ing that one spin relaxation mechanism dominates only. 
Therewith in spite of the strong anisotropy of spin-orbit 
scattering, the relaxation times of spin lying in the plane 
of a heterostructure with zinc-blende lattice turn out to 
be independent on orientation with respect to crystallo- 
graphic axes. 

Real heterostructures differ from investigated ideal 
objects in that several spin relaxation mechanisms ex- 
ist [^|4|. A spin relaxation mechanism due to only cubic 
in wave vector terms of the bulk Hamiltonian was inves- 
tigated for rectangular quantum wells (QWs) in Ref. [g. 
It was noted that even in an asymmetrical GaAs QW, the 
efficiency of another mechanism due to linear in two di- 
mensional wave vector terms is negligibly small. The au- 
thors of Ref. |0| analysing experimental data on anoma- 
lous magnetoresistance in InGaAs QWs demonstrated 
that the both mechanisms may be comparable in effi- 
ciency. But in Refs. ^j^, it was mentioned that the both 
mechanisms are additive in spin relaxation. 

This communication is devoted to an investigation of 
spin relaxation processes when several mechanisms of 
spin-orbit scattering exist. We show that contributions 
of these mechanisms interfere and their simultaneous ac- 
tion leads to the strong anisotropy of spin relaxation even 
in the plane of a QW. 



II. THEORY 

In zinc-blende semiconductors, spin relaxation of elec- 
trons is well known to be due to spin-orbit splitting of 
conduction band. In a bulk crystal, the splitting is cubic 
in wave vector. In a QW structure, the corresponding 
Hamiltonian has to be averaged over the motion along 



the growth axis. We consider the QW grown along z- 
direction parallel to [001] and choose x and y directions 
coinciding with crystallographic axes. At relatively small 
carrier concentrations, one can neglect cubic in 2D wave 
vector terms and the Hamiltonian has the form: 



Hi = ai{axkx - CTyk, 



yrvyj 



(1) 



Here ai {i = x, y) is the Pauli matrix, ki is the wave 
vector component in the plane of the QW and ai is a 
constant. 

In asymmetrical heterostructures, there is a contribu- 
tion to the Hamiltonian which is absent in the bulk [o|: 



H2 = a2{axky - Uykx) , 



(2) 



where 02 is the constant determined by heterointerface 
properties. 

To calculate spin relaxation times, one has to consider 
the dynamics of spin density matrix, p: 



where the total two dimensional Hamiltonian is: 



H 



2m 



V + H' . 



(3) 



(4) 



Here m is an effective electron mass, V{x, y) is a scatter- 
ing potential and 

H' = Hi + H2. 

We assume that the scattering is elastic and independent 
on spin indices. 

Since H' is a small perturbation, the spin relaxation 
times turn out to be much longer than isotropisation 
times of momentum distribution of electrons. For this 
reason it is convenient to represent the density matrix as 
a sum W, : 



p^ p + p 



0. 



where the bar means averaging over the directions of k. 



Here p depends on e — h k /2m and the anisotropic 
part of the density matrix is due to H' only. Hence p' is 
proportional to H' , i.e. p' is small in coinparison to p. 
Then in the first order in H' , the Eq. (|^) has the form: 



%--^ [H',p] J2 ^kk' [p'ik] - p'(fc')] 



l-s^ 



(5) 



(6) 



Here Wkk' is the probability for scattering from the po- 
tential V from the state with k to the state with k' . 

One can see from the Eq. (g) that p' relaxes in the 
time which is of order of momentum relaxation time, but 
p relaxes in the longer time which is determined by H' 
— see (^). At these long times, the Eq. (||) reduces to a 
quasistationary one: 



Y,W^^,[p'{k)~p'{k')]^~^[H',p] 



h 



(7) 



Finding p'{k) from here and substituting it to the 
Eq. (p), one can obtain the closed equation for p: 



i = -^E-"[^-«.K.p]] 



(8) 



Here 



K^j'^ H'{k) cxp(-zn(pk) , (9) 

where tpk is the angle between k and x axis, and 

— ^(bde VFkk'(l - cos ri6l) , (10) 

T« J 

where — (fi^ — (pk' ■ The Eq. (^) clearly demonstrates 
that it is spin-orbit interaction which causes p relaxation. 
After substituting p in a form 

P = /o + 2*^ ■ ^ ' 
the Eq. (^) reduces to following equations: 

foie,t)=0, (11) 



sei{e,t) = -Tij{e)iSj{e,t) 



where 



^=-^Tr|^T„[i/i„,K,a,]]a,l 



(12) 



(13) 



An initial condition may be derived considering the 
spin dynamics after the time t„, but before the spin re- 
laxation time. In the time t„, the density matrix becomes 
isotropic but the spin relaxation processes do not start 
yet. Therefore |]: 



/o(e) = -[F+(e)-HF_(e)] 



ae(e) = s[F+(£)-F_(e) 



(14) 



(15) 



where s is the unit vector along the spin and F±{e) are 
distribution functions of particles with the spin projec- 
tion on s equal to ±1/2. 



Taking into account that the spin density, S{t), is the 
average of ae over s, one can obtain the kinetic equation 
for it at the time longer than t„: 



S, = -^ 



(16) 



where the tensor of reciprocal spin relaxation times is: 
1 _Jde[F+{e)-F_{e)]r,,{e) 



Jde[F+{e)~F^ie)] 



(17) 



The Eq. (|l^) represents the extension of the results of 
Ref. [pj to the case of an arbitrary spin-orbit interaction 
H' and takes into account the anisotropy of scattering. 

Substituting H' = Hi + H2 into (|l^ ) and then Tij 
into (p7|) we have: 



±^C{al + al),— = — = Q., 



(18) 



lie 1 

= = 17 (a? +02), = -Caia2, (19) 



' XX > yy 
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where 



C 



8m / de [F+{e) - F_ (e)]Ti (e) e 
^ Jde[F+{e)-F_{e)] 



(20) 



Equations (|l^ — ^) generalize the results ^,^,|| for 
the case of two spin relaxation mechanisms. 



III. DISCUSSION 

It follows from the Eq. (nm that because of the two 
mechanisms (oi -02 7^ 0), the spin relaxation in the plane 
of the QW becomes anisotropic. It should be empha- 
sized that if there is only one mechanism (ai -02 — 0), 
then the spin relaxation is isotropic in spite of the cubic 
symmetry of the Hamiltonian Hi or H2 ■ Thus the cubic 
anisotropy of the real QW structure manifests itself due 
to the interference of two spin relaxation mechanisms. 

The system (|l^) may be rewritten as follows: 



where 



Q I Q _ _^x^^y 

^x ^ ^v — 



1 C-/ n2 

— = 17(01 ±02) 

T± 2 



(21) 



(22) 



The times t+ and r_ describe the relaxation of the spin 
oriented along the directions [110] and [110] respectively. 
The most bright manifestation of spin relaxation 
anisotropy occurs if ai = ±a2. In this case, one of the 
times r-|_ or r_ becomes infinite. Therewith the other is 
equal to t^^. 



Efficiency of the mecfianisms depends on both the ma- 
terial and the shape of the QW. It was shown that in a 
rectangular GaAs/AlGaAs QW, the mechanism (|]) dom- 
inates [g[, and in an asymmetrical InGaAs/AlAs QWs 
the mechanism (0) is the most important [nO| or they are 
comparable Q. 

It is clear that if both ai and 02 are not equal to zero, 
then the spin sublevels split. In the work ||ll[| it is shown 
that oi and 02 may be comparable in magnitude and, 
hence, the spin splitting is strongly anisotropic. 

The spin relaxation anisotropy results from the initial 
Td symmetry of the zinc-blende semiconductor. For this 
reason, the similar effect can take place in a strained 
bulk crystal. The corresponding Hamiltonian linear in 
3D wave vector, k, and components of an elastic strain 
tensor, Ui,, has the form: 



H'{u) = Ai Uii (tJj+ifci+i - ai+-2.ki+2) 
+A2 Uij {aikj - CFjki) . 



(23) 



Here i, j — x,y, z, i + 3 ^ i, Ai and A2 are constants. 
Doing calculations for 1/tjj in a way similar to (p^), one 
can obtain three different spin relaxation times. It can 
be shown that the maximum anisotropy may be achieved 
if 



AlUxx = MUyy = -AlUzzl'i. = A2Uxy/'i 



(24) 



with the rest of Uij — 0. Therewith two spin relaxation 
times are equal to each other and the third is infinite. 
Note that the tensor Uy determined by (p3) may be ob- 
tained by applying two uni-axial strains along the axes 
[001] and [110] and they are not restricted to uni-axial 
strain along any axes. 
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IV. CONCLUSION 

The possibility for spin relaxation suppression was 
noted in Ref. p] for a QW grown along [110] direction 
when the spin is oriented along the same axis. The 
present work shows that the spin relaxation rate also de- 
creases for [110] direction, but in a QW grown in the sym- 
metrical direction [001]. Therefore this decrease takes 
place for the spin lying in the plane of the QW. 

Analysing weak localisation effect, the authors of 
Refs. [HJj] showed that the mechanisms (|l|) and (||) sup- 
press each other in anomalous magnetoresistance, but 
they are additive in spin relaxation. The present analy- 
sis shows that the suppression occurs in the spin relax- 
ation also. Besides, we have found that spin relaxation 
is anisotropic even in the plane of the QW. 
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